Width of the — rr phase transition in diffusive magnetic Josephson junctions 
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We investigate the Josephson current between two superconductors (S) which are connected 
through a diffusive magnetic junction with a complex structure (F c ). Using the quantum circuit 
theory, we obtain the phase diagram of and n Josephson couplings for F c being a IFI (insulator- 
ferromagnet-insulator) double barrier junction or a IFNFI structure (where N indicates a normal 
metal layer). Compared to a simple SFS structure, we find that the width of the transition, defined 
by the interval of exchange fields in which a — n transition is possible, is increased by insulating 
barriers at the interfaces and also by the presence of the additional N layer. The widest transition is 
found for symmetric F c structures. The symmetric SIFNFIS presents the most favorable condition 
to detect the temperature induced — tt transition with a relative width, which is five times larger 
than that of the corresponding simple SFS structure. 

PACS numbers: 74.45.+C, 74.50.+r, 72.25.-b 
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I. INTRODUCTION 

Ferromagnet-superconductor (FS) heterostructures 
feature novel and interesting phenomena, which have 
been active topics of investigation for more than half a 
centurjsii2^. Meanwhile, Josephson structures compris- 
ing a ferromagnetic weak link have been studied exten- 
sively. The existence of the 7r-junction in a SFS struc- 
ture is one of the most interesting phenomena which oc- 
curs for certain thicknesses and exchange fields of the 
F layer 7 - 8 - 9 - 10 - 11 - 12 - 13 - 14 - 15,16 - 17,18 - 19 - 20 - 21 - 22 - 23 This man- 
ifestation is due to the oscillatory behavior of the su- 
perconducting pair amplitude and the electronic density 
of states in the ferromagnei^^. In a 7r-junction the 
ground state phase difference between two coupled su- 
perconductors is 7T instead of as in the usual 0-state 
SNS junctions. The existence of a 7r-state was predicted 
theoretically by Bulaevski et alZ and Buzdin et al£, and 
has been first observed by Ryazanov et al£. The tran- 
sition from the 0- to 7r-state is associated with a sign 
change of the critical current, I c , which leads to a cusp- 
like dependence of the absolute values of I c on tempera- 
ture. Later, the nonmonotonic temperature dependence 
of the critical current in diffusive contacts was observed 
in other experiment o 10 ! 11 ' 12 ! 13 ' and was attributed to 
the — 7T transition induced by the ferromagnetic ex- 
change field. The — tt transition has been studied 
theoretically by several authors in clea n 15 i 16 i 17 i 18 and 
diffusiv o 17 i 18 i 19 i 20 i 21 i 22 i 23 Josephson contacts for different 
conditions and barriers at the FS interfaces. 

An interesting application of a 7r-junction is a super- 
conducting qubit as one of the most noticeable candidates 
for a basic element of quantum computing. Further- 
more, 7r-junctions have been proposed as phase qubit ele- 
ments in superconducting logic circuit o 24 ' 25 ! 26 ' 27 ! 28 . Also, 
a phase qubit in SIFIS junctions, in which the qubit state 
is characterized by the and the tt phase states of the 
junction, has recently been suggested^ 9 .. Due to these 
exciting proposed applications, the detection of — tt 



transitions with very high sensitivity is necessary. Inves- 
tigating the phase diagrams of — tt transitions 1 ^ , for 
different structures with different characteristics should 
make it possible to determine the most efficient control 
of the — tt transition. 

In this paper we investigate the width of the 
temperature-induced — tt transition in a diffusive SF C S 
junction. Here, F c represents a complex ferromagnetic 
junction of length L, which consists of diffusive ferromag- 
netic and normal metallic parts as well as insulating bar- 
riers. We define the width of the transition as the interval 
of exchange fields, in which the temperature-dependent 
transition from the 0- to 7r-phase is possible. We use 
the so-called quantum circuit theory (QCT), which is a 
finite-element technique for quasiclassical Green's func- 
tions in the diffusive limit^^ 2 .^. The QCT-description 
has the advantage, that it does not require to specify 
a concrete geometry. By a discretization of the Usadel 
equation 3 ^ one obtains relations analogous to the Kirch- 
hoff laws for classical electric circuit theory. These re- 
lations can be solved numerically by iterative methods 
and one obtains the quasiclassical Green's function of the 
whole system. The QCT has been generalized to spin- 
dependent transport in Ref. [35ll36l . We adopt the find- 
ing of this paper for FN contacts to handle our problem 
of the SF C S contacts. We discretize the interlayer be- 
tween the superconducting reservoirs into nodes. Follow- 
ing Refs. 36.37:, every node in a ferromagnetic layer with 
specific exchange field is equivalent to a normal node con- 
nected to a ferromagnetic insulator reservoir which deter- 
mines the exchange field. This similarity has been ver- 
ified experimentally with EuO|Al|Al203|Al junctions 38 . 
It has been found that the induced exchange-field of the 
EuO-insulator, which is responsible for spin-splitting in 
the measured density of states, was of the same order 
as its magnetization^. Also, the authors of Ref. have 
shown that the normalized density of states in the normal 
metal, which is connected to a superconductor and an in- 
sulator ferromagnet at its ends, is the same as the one for 
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a BCS superconductor in the presence of a spin-splitting 
magnetic fiel d 39 ' 40 . This method allows us to calculate 
the Josephson current flowing through the SF C S contact 
for an arbitrary length L and all temperatures while fully 
taking into account the nonlinear effects of the induced 
superconducting correlations. 

We investigate the width of the transition, Ah, for four 
different cases of SF C S structures with ideally transpar- 
ent FS-interfaces, symmetric SIFIS and asymmetric S1FS 
double barrier F-junctions, and more complicated SIFN- 
FIS structures (where I and N denote, respectively, insu- 
lating barrier and normal metal). For a fixed length L, all 
these systems show several transition lines in the phase 
diagram of T/T c and h/T c . Higher order transitions oc- 
cur at large exchange fields h. We find that higher order 
transitions are wider than the first transition. Also, de- 
creasing the contact length L leads to a widening of the 
transitions and, at the same time, to an increasing of the 
exchange field, hi n , at which the transition starts. Nev- 
ertheless, the relative width of the transition, given by 
the ratio Ah/hi n , decreases. 

For the SIFIS structure we show that the existence of 
the I-barriers at the FS interfaces broadens the — it 
transitions and, hence, improves the conditions to detect 
such transitions. In addition, we find that a symmet- 
ric double barrier structure with the two barriers having 
the same conductance shows wider transitions than the 
corresponding asymmetric structure with the same total 
conductance but different conductances of the barriers. 
An even larger width of transitions can be achieved by 
including an additional normal metal part into F c . This 
motivates our study of an SIFNFIS structure, for which 
relative width Ah/hi n is in general larger than that of 
the corresponding SIFIS with the same total conductance 
and the mean exchange field of the F c part. 

The structure of this paper is as follows. In Sec. [Til 
we introduce the model and basic equations, which are 
used to investigate the SF C S- Josephson junction. We in- 
troduce the finite-element description of our structures 
using quantum circuit theory technique. In Sec. Hill wc 
investigate phase diagrams of — it transitions for the 
SFS, SIFIS, and SIFNFIS structures. Analyzing our find- 
ings, we determine the most favorable conditions for an 
experimental detection of the — it transitions. Finally, 
we conclude in Sec. [TV] 



II. MODEL AND BASIC EQUATIONS 

We consider a ferromagnetic SF C S Josephson structure 
in which two conventional superconducting reservoirs are 
connected by a complex diffusive F c junction. We in- 
vestigate three cases of F c : (i) a simple F layer with a 
homogeneous spin-splitting exchange field h (SFS), (ii) 
a double barrier IFI structure, in which the F-layer is 
connected via I-barriers to the reservoirs (SIFIS), and 
(iii) an IFNFI junction composed of two ferromagnetic 
layers with the same length Lp and the same exchange 
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FIG. 1: (Color online) Quantum circuit model of the SF C S 
structure. The contact F c is discretized into n nodes which 
are connected to each other by tunnel barriers of conductance 
gr\ gsiF and qfs2 denote the conductances of S1F and FS2 
interfaces, respectively. The inverse of the average level spac- 
ing, 5, represents the leakage term due to a finite volume of a 
node; FIR represents a ferromagnetic insulating reservoir. 

field and a normal metal with length Ln in between such 
that L = L N + 2L F (SIFNFIS). We compare the width of 
the temperature- induced — tt transitions for these three 
types of structures. In all cases F c has the same length, 
total conductance, and the mean exchange field h. 

In our approach, we make use of the quantum cir- 
cuit theory which is a finite-clement theory technique 
for quasiclassical Green's function method in diffusive 
limi t 31 ' 32 ' 33 . In this technique, each part of the struc- 
ture is represented by a node which is connected to 
other nodes or superconductor /ferromagnet reservoirs 3 ^. 
Green's functions are calculated by using balance equa- 
tions for matrix currents entering from the connectors, 
which is described in terms of its transmission properties 
and Green's functions of the nodes forming it, to each 
node. For calculations we follow the procedure similar to 
that of Ref. |37|. We discretize the conducting part of F c 
into n nodes as presented in Fig. [1] A node in the ferro- 
magnetic part will be presented by a normal-metal node 
connected to a ferromagnetic insulating reservoir (FIR) 
which induces an exchange field equal to the exchange 
field of the ferromagnetic part at the place of the node. 

Each of the superconducting reservoirs is assumed to 
be a standard BCS-superconductor. Our circuit connect- 
ing those reservoirs consists of different types of nodes 
in F c . One type are the normal nodes in the middle 
of F c , each of which is connected to two neighboring 
nodes which are either normal nodes or F nodes. An- 
other type are F nodes placed at the two ends of F c , 
where each of them, in addition to its connection to 
two neighboring nodes, is connected to FIR as well and, 
hence, feels the exchange field directly. As can be seen 
in Fig. [TJ each of the two neighboring nodes of a F node 
can be another F node, an N node, or a superconducting 
node. We denote the conductances of the tunnel bar- 
riers at S1F and FS2 interfaces by gsiF and gFS2, re- 
spectively. Also, gr represents the conductance of the 
tunnel barrier between each two nodes inside F c ; gx is 
determined by gp a , the total conductance of F c , exclud- 
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ing the conductances of the barriers at the interfaces 
(n - l)/g T = (1/.9fJ - (I/ssif + l/fl^sra)- In general, a 
node i is characterized by a Green's function Gi, which 
is an energy-dependent 4 x 4-matrix in the Nambu and 
spin spaces. Furthermore, all nodes in F c are assumed to 
be coupled to each other by means of tunneling contacts. 
However, a finite volume of a node and the associated 
decoherence between electron and hole excitations are 
taken into account by the leakage matrix current which 
is proportional to the energy, e, and the inverse of the 
average level spacing in the node, < 5 31 ' 32 . 

For a structure with spin-dependent magnetic contacts 
and in the presence of F and S reservoirs, the matrix cur- 
rent was developed in Ref. [36|. In the limit of tunneling 
contacts, which is our interest, the matrix current be- 
tween two nodes i,j is defined a o 36 i 41 , 
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The first term demonstrates the usual boundary condi- 
tion for a tunneling junction, where gij is the tunneling 
conductance of the contact between the two nodes. The 
second term exists due to the different conductances for 
different spin directions which leads to the spin polar- 
ized current through the contact. We assume this term 
to be negligible as, G M r ~ gjj - g\j < 9i,j- Also, 

Gq = e 2 /2irh is the quantum of conductance, h is the 
exchange field of the node, and a and r are the vectors 
consisting of Pauli matrices in spin and Nambu space. 

Using Eq. ([I]) for different matrix currents entering 
into a given node i, we apply the condition of current 
conservation to obtain the following balance equation, 



[ E 

j=i— 



9j, 



.G 



G c 



Gj + i—(hi.a)f 3 - i—eT 3 a ,Gi} = 0. 



(2) 

Here, the first term represents the matrix currents from 
neighboring nodes i — l,i + 1, which could be F, N or 
S. The second and third terms are, respectively, the ex- 
change term and the leakage matrix current. Also, <jq 
represents unit matrix in spin space. 

We consider the case, in which the exchange field in the 
ferromagnetic parts of F c is homogeneous and collinear. 
Then, it is sufficient to work with the 2x2 matrix Green's 
function of spin-er (a =f / j) electrons in Nambu space. 
In the Matsubara formalism the energy e is replaced by 
Matsubara frequency ioj — iTrT(2m+l)) and the Green's 
function has the form 



G = 



G F 

F* -G 



(3) 



Neglecting the inverse proximity effect in the right and 
left superconducting reservoirs, we set the boundary con- 
ditions at the corresponding nodes SI and S2 to the bulk 



values of the matrix Green's functions: 



Vuj 2 + A 2 V Ae^/ 2 



(4) 



Here Ae ±J */ 2 are, respectively, the superconducting or- 
der parameters in the right and left superconductors, 
and 4> is the phase difference. The matrix Greens func- 
tion satisfies the normalization condition, G — 1. The 
temperature-dependence of the superconducting gap A 
is modeled by the following formula 4 —^ 3 - 



A = !.7(ir,.tanli( l.71\/— - I). 



(5) 



We scale the size of F c in units of the diffusive su- 
perconducting coherence length, £5 = \J^li m p where 
Co = vp/irAo with vf being the Fermi velocity and 
Aq = A(T = 0) = 1.76T C , and li mp is the mean free 
path in the F-layer related to the diffusion coefficient 

via D — v F kmp/3. Two more scales that we use 
are h/T c and T/T Cl where T c is the critical tempera- 
ture of S reservoirs. Also, the mean level spacing de- 
pends on the size of the system via the Thouless energy 
Etk = D/L 2 = gTS/(n — 1)Gq (Planck and Boltzmann 
constants, Ti and Icb, are set to 1 throughout this paper). 

In the absence of spin- flip scatterings, the balance 
equation, Eq. ([2]), is written for each spin direction sep- 
arately for all n nodes in F c . This results in a set of 
equations for n matrix Green's functions of the nodes 
that are solved numerically by iteration. In our calcu- 
lation we start with choosing a trial form of the matrix 
Green's functions of the nodes, for a given cf>, T, and the 
Matsubara frequency m = 1. Then, using Eq. ^ and 
the boundary conditions iteratively, we refine the initial 
values until the Green's functions are calculated in each 
of n nodes with the desired accuracy. Note that in gen- 
eral for any phase differennce 0, the resulting Green's 
functions vary from one node to another, simulating the 
spatial variation along the F c contact. From the result- 
ing Green's functions we calculate the spectral current 
using Eq. fTJ and obtain 



T 



U) m —00 
LU m — — 00 



(6) 



In the second step we set the next Matsubara frequency 
m = 2, find its contribution to the spectral current, and 
continue to the higher frequencies until the required pre- 
cision of the summation over m is achieved. Finding the 
spectral current, for the given temperature and phase 
difference, enables us to obtain the dependence of the 
critical current I c on T. Finally, we increase the number 
of nodes, n, and repeat the above procedure until all the 
spectral currents for every temperature and phase differ- 
ence reach the specified accuracy. We find that for typical 
values of the involved parameters, a mesh of 60 nodes is 
sufficient to obtain I c through the diffusive F c structure 
with an accuracy of I0 -3 across the whole temperature 
range. 
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III. RESULTS AND DISCUSSIONS 

/From the numerical calculations, described above, we 
have obtained the phase diagram of 0— tt transition in the 
plane of h/T c and T/T c . We analyze the width of — tt 
transitions for the SFS, symmetric SIFIS and asymmetric 
SIFS double barrier junctions, and SIFNFIS structures. 

Concerning such transitions, the width Ah defines the 
interval of the h, in which a temperature-induced tran- 
sition occurs. We compare relative width, the ratio 
Ah/hi n , of different structures, where hi n is the exchange 
field in which the transition starts (see Fig. [2^). In prac- 
tice, we fix the size of the structures, L/£ s , and then vary 
the value of h/T c for detecting the change in the sign of 
the critical supercurrent as the transition occurs. We ex- 
pect that the detection of a — tt transition can be more 
feasible for the structure having larger Ah/hi n . 

A. SFS structures 

First, we consider the SFS structure. Figure. [2^ 
presents the typical 0-7r transitions for such a junction 
with L/£ s = 1.5, where the supercurrent is scaled in units 
of Io = (7r/2)Ao/ei?F e - Here, Rp c is the total resistance 
of F c . We observe that the nonzero supercurrent at the 
transition point is larger when the transition temperature 
is lower. Also, the phase diagram is shown in Fig. ^Bp in 
the vicinity of the first and the second — tt transitions. 
At the first transition the junction goes from the 0- to 
the 7r-state starting at hi„ and T = 0. Increasing h, the 
transition temperature increases toward T c and above the 
value h = hi n + Ah, the junction will be in tt state at 
all temperatures. Increasing h further, the junction stays 
at its tt state until the exchange field reaches the value 
at which the second transition starts (see Fig. [2Jd) where 
the junction changes back to a 0-state. In principal, it is 
possible to go to the higher exchange fields to see higher 
transitions. However the amplitude of the supercurrent 
will be extremely small and difficult to detect experimen- 
tally. 

We have observed that the second transition is always 
wider than the first one. In the case of Fig.[5j3, the width 
of the first transition is nearly 0.65 of that of the second 
one. Furthermore, the relative width for first transition 
is 0.20, while the second transition has Ah/hi n = 0.06. 
This finding can also be generalized to higher transitions. 
In brief, higher transitions are always associated with 
larger widths. In spite of having a smaller width, the 
first transition seems to be more feasibly detectable, since 
they have higher Ah/hi n . 

Looking at the origin of the existence of — tt tran- 
sition, we can understand this finding. An oscillating 
behavior of the order parameter in a ferromagnetic layer 
makes the occurrence of different signs of order parame- 
ters of the superconductor reservoirs, possible. This ef- 
fect, being in charge of the 7r-phase state, can be seen 
when the length of the ferromagnet is of the order of 
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FIG. 2: (a) (Color online) Normalized critical current, I c /Io, 
versus temperature for a SFS structure with the length 
Lf/£s = 1.5 and different exchange fields h/T c . The inset 
shows the current-phase characteristic for h/T c = 5.1 in the 
vicinity of the — n transition temperature, (b) The corre- 
sponding — tt phase diagram showing the phase boundaries 
up to the second transition, (c) Phase diagram of the first 
— tv transition for different length of the junction Lf/£ s . 



half-integer multiple of a period 2tt£f, where £p is the 
ferromagnetic coherence length of the ferromagnet. In 
the diffusive limit when h >> T c , this length is equal to 
£f = \J D jh. Hence, as d^p/dh is inversely proportional 
to the exchange field, when the exchange field becomes 
larger the rate of reduction of £f decreases and the sys- 
tem will remain longer in the region of transition. 

Now let us consider the effect of the length of F on 
the width of the transition. In Fig. [2J;, the width of 
the first transition for three lengths L/^ s = 0.5,1,2 are 
compared. As mentioned above, the condition for the 
occurrence of the first transition is that the length L 
becomes of the order of half-integer of the period. For 
a smaller L this condition is fulfilled at larger h which, 
in light of the above discussion, means a wider — tt 
transition. This can be seen easily in Fig. [2fc. Note that 
the transition between the two states always starts from 
lower temperatures. 
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FIG. 3: (a) (Color online) Normalized critical current, 
Ic/Io, versus temperature of a symmetric SIFIS structure with 
gsiF = 9FS2 = 0.018(?t for different h/T c , when Lf/£, s = 1-5. 
The inset shows the corresponding I — <f> characteristic for 
h/T c = 2.0 in the vicinity of the — n transition temper- 
ature, (b) Phase diagram of — tt transition for a sym- 
metric SIFIS structure with gsiF = gFS2 ~ 0.018<7t (solid 
line) and the corresponding asymmetric structure (dashed 
line) with gsiF = 0.1(?fs2 = O.lgr, when Lf/£s = 1.5. (c) 
The same as (b) but for a symmetric SIFNFIS structure of 
gsiF = 9FS2 = 0.018<7t with L/£ s = 1-5 and various Ljv/Lf- 



B. SIFIS, SIFS, and SIFNFIS structures 

Next, we examine the effect of putting insulating bar- 
riers at FS-interfaces. In Fig. the typical 0-7T transi- 
tions for SIFIS structure with L/£ 3 = 1.5 is shown. As 
one can see, the presence of barriers adjusts the nonzero 
minimum cusp appearing in the diagram of the criti- 
cal current versus temperature. Figure. [3Jd manifests 
the — 7r phase diagram for the symmetric SIFIS (solid 
curve) and the asymmetric SIFS (dashed curve) double 
barrier F c . Compared to the corresponding SFS with 
Ah/hin — 0.2, these structures show wider transitions 
with Ah/hi n — 0-93 for SIFIS of the conductance of the 
barrier gsiF = 9FS2 — 0.018gT, and Ah/hi n = 0.61 for 
SIFS of g S iF = 0.1gFS2 = O.lgT- 

We have found that the strength of the barriers be- 
tween the FS junctions is the most important parameter 



for determining the width of 0— tt transitions. On the one 
hand, as the barriers get stronger the width of transitions 
becomes wider. This widening will be more pronounced 
for short length structures. On the other hand, for these 
structure the transition will start from a lower exchange 
field in comparison with the corresponding SFS systems. 

Considering the effect of the relative values of the con- 
ductances of the two barrier, a symmetric SIFIS struc- 
tures shows broader transitions as compared to the asym- 
metric SIFS structure with the same total conductance, 
as can be seen in Fig. [5b. 

In addition, considering the displacement of the bar- 
rier in a SlIiFiI 2 F 2 S2 hybrid structures, we have found 
that the effect of barriers becomes more important as 
the barriers are closer to the ends of F c , so that, SIFIS 
is the most optimal structure regarding the width of the 
transitions. 

Finally, we have investigated the width of the — tt 
transitions for SIFNFIS structures. The phase diagrams 
are shown in Fig. [3J: for junctions with L/£ s = 1.5 and 
various values of the length of the N part, Ljy. We see 
that, putting a normal metal between the ferromagnets 
while keeping the magnetization of the system constant 
increases the width of the transition somewhat. This can 
be due to stronger penetration of superconductivity near 
the FS boundaries where the density of magnetization 
is larger which strengthens the mean effect of exchange 
field. 

We have also observed that increasing Lm leads to a 
further increase in the width of transition. However, this 
increase is saturated at higher lengthes. While the width 
for SIFNFIS structures of Lm = 2Lp is almost doubled 
compared to the SIFS structure, it is increased only few 
percent by increasing from 2Lp to 4Lp (see Fig.|3J;). 

It is worth to note that taking the absolute width Ah 
as measure of the feasibility to detect the temperature- 
induced — tt transition will lead to similar results as 
those of obtained above by considering the relative width 
Ah/hi n . However, the definition by Ah/hm seems to be 
more appropriate, since higher feasibility of detection re- 
quires not only larger Ah, but also smaller hi n in order 
to have weaker exchange-induced suppression of the su- 
percurrent. 



IV. CONCLUSION 

In conclusion, we have investigated the width of — tt 
transitions for various diffusive ferromagnetic Josephson 
structures (F c ) made of feromagnetic (F) and normal 
metal (N) layers and the insulating barrier (I) contacts. 
We have performed numerical calculations of the Joseph- 
son current within the quantum circuit theory technique 
which takes into account fully nonlinear proximity effect. 
The resulting phase diagram of and tt Josephson cou- 
plings in the plane of T/T c and h/T c shows that the ex- 
istence of the insulating barrier contacts and the normal 
metal inter-layer leads to the enhancement of the relative 
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width of the temperature induced transition. The rela- 
tive width is parameterized by the ratio Ah/h in with Ah 
and hi n being, respectively, the exchange field interval 
upon which the transition is possible and the initial value 
of h at which the transition occurs at T = 0. We have 
also observed that while the conductance, the magneti- 
zation, and the length of the F c junction are kept fixed, 
symmetric structures with the same barrier contacts and 
the same F layers in a SIFNFIS structure show larger 
relative width of the transition compared to that of the 
asymmetric structures. Among the studied structures, a 
symmetric SIFNFIS junction have the highest Ah/hi n , 



which makes it more practicable for highly sensitive de- 
tection of the temperature-induced — it transition. 
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